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In high energy hadron-hadron and e + e~ collisions, to isolate a part of the phase space in multi- 
hadron final states is necessary for exploring the underlying dynamics. It is shown that the partition 
of phase space according to the value of rapidity, popularly used in hadron-hadron collisions, is 
inappropriate for the study of e + e~ collisions. The proper way in the latter case is to identify the 
visible jets and take them as objects for detailed study, forming an extended phase space. The value 
of the distance-measure for the identification of visible jets is determined. A new variable r is 
introduced to further partition the phase space inside jets, which possesses a very good anomalous 
scaling property, showing that the r-distribution inside jets is a self-similar fractal. 
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The dynamics of multi-hadron final states in high en- 
ergy collisions is an idea laboratory for the study of the 
perturbative and nonperturbative properties of the basic 
theory of strong interaction — QuantumChromo Dynamics 
QCD. This study, starting from the sixties of last cen- 
tury, has a long history. The multiplicity as well as the 
phase space volume have increased for orders of magni- 
tude since then. It becomes necessary to take a part of 
the phase space for detailed analysis instead of to merely 
study it as a whole. 

In hadron-hadron collisions, already in the 1970s the 
ISR experiments started to study the "central rapidity 
region" \y\ < Y cut 0, where the rapidity y is defined as 



y = - In — 
y 2 E-p 



P\\ 



(1) 



with E the energy and p\\ the longitudinal momentum of 
the particle in the cm. frame, Y^. ut is a specially chosen 
cut smaller than the maximum rapidity reached in the 
collisions. The longitudinal direction is taken to be the 
direction of the momenta of the two incident particles. 

Later, the successful operation of the p-p collider SppS 
at CERN raised the cm. energy </s for an order of mag- 
nitude. It enabled the UA5 Q and other collaborations 
to carry on their study in various rapidity windows. 

More recently, this approach has been extended to 
e + e~ collisions. For example, Delphi || and opal || 
Collaborations have studied the intermittency and cor- 
relations in hadronic Z° decay in a restricted rapidity 
region \y\ < 2, with the thrust or sphericity axis taken as 
the longitudinal direction. 

However, the dynamics of c + e~ collisions is different 
from that of hadron-hadron collisions. In hadron-hadron 
collisions the central rapidity region has a special phys- 
ical meaning. Particle production occurs mainly in this 
region. The regions outside of the central one are the 
fragmentation regions, where the leading particle effect is 
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present. On the other hand, in e + e~ collisions the final- 
state system comes from a point source (virtual photon 
or Z°). There is no leading particle effect, or, if some 
similar effect exists in the individual jets, it is not di- 
rected along the thrust (or sphericity) axis in multi-jet 
events, see e.g. Fig.l and Fig. 2. In this case, to study 
the "central rapidity region" with thrust or (sphericity) 
axis as longitudinal direction is doubtful. 

In 1994 J. D. Bjorkcn [|| proposed to isolate the gluon 
jets for the study of the crucial features of the underlying 
dynamics and put forward the idea of "extended phase- 
space", cf. figure 1 of Ref. In this approach gluon jets 
produced in e + e~ collisions are first identified and then 
taken out as a part of the extended phase-space. Further 
partition of phase space is done inside the jets. 

Different ways of cutting phase-space, as e.g. those 
done in Ref . || || and those proposed in Ref . (5| , will ap- 
parently result in different physics. What is the proper 
way for the partition of phase-space? How to realize 
it? are urgent problems not only for the study of e + e _ 
collisions but also for the study of the ultra-high en- 
ergy hadron-hadron and/or relativistic heavy ion colli- 
sions at Fermilab- TE VAT RON, Brookhaven-RHIC and CERN- 
lhc, where a vast amount of jets are being or will be 
produced. 

In this letter we will try to give an answer to these ques- 
tions. For illustration, we will use LUND Monte Carlo JET- 
SET?^ to generate a sample of 1.5 million c + c~ collision 
events at ^/s = 91.2 GeV. 

Let us first analyse the usefulness of "central rapidity 
region" in e + e~ collisions. As already pointed out, the 
final state hadrons in these collisions are produced from 
a point source. There is no leading particle effect, or in 
multi-jet events such effect is not along the "longitudi- 
nal" (thrust or sphericity) direction. In these collisions 
there are 2-jet, 3-jet, 4-jet, ... events. A cut in rapid- 
ity (with the thrust or sphericity axis as "longitudinal") 
will have very different effect on different-number-of-jet 
events. For example, in a 2-jet event a y cut will cut out 
the most energetic particles in both jets; while in a 3-jet 
event a y cut will only cut out the most energetic parti- 
cles of one jet but retain almost all the particles of the 
other two jets. Two typical examples are given in Fig.l, 
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FIG. 1: The lego plot of two typical 3-jet events determined 
by DURHAM algorithm with log d cut = —2.2. Downward and 
upward triangles are the particles of the two right-hand-side 
jets, respectively. Open stars are the corresponding jet axes. 
Dashed circles have radius r = 0.7. 







where the dotted lines denote the rapidity cut \y\ < 2. It 
cuts out a large part of the left-hand-side jet but retains 
almost all the particles of the other two jets. The results 
from such an asymmetric cut is hard to be interpreted. 

Physically, a more appropriate approach is not to se- 
lect a group of particles with rapidity in some range but 
to isolate the individual jets. The corresponding cut con- 
dition is then r < R cu t |§ instead of y < Y cut , where Q 



\J iv - Vict) 2 + (<P- ^jot) 5 



(2) 



rj and (p are the pseudo-rapidity and azimuthal angle of a 
final state particle in a jet, respectively, while ?7j c t and <pj C t 
are those of the jet. i? C ut is a number around unity. For 
the dashed circles in Fig.l this number has been chosen 
to be R cut = 0.7 §. Note that in Eq.(2) we have used 
pseudo-rapidity rj defined as 



■q = -lntan(0/2) 



(3) 



instead of the rapidity defined in Eq.(l). At high energy, 
■q w y. For the study of phase space partition, rj-ip is 
a natural choise, which are directly connected with the 
corresponding variables in a space-time description ||. 
Note also that, when a jet has its axis (the total momen- 
tum of the particles inside it) nearly collinear with the 
thrust (or sphericity) axis, which is used as longitudinal 
to define the variables -q and ip, the azimuthal angle ip 
of the particles inside the jet will be spreading and the 
jet will not be visible in the rj-tp lego-plot, cf. the left- 
hand-side jets of Fig.l and Fig's. 2 (a), (c). Therefore, 
in the following we will restrict ourself to consider only 
those jets whose axes have an angle greater than 15° and 
less than 165° with the thrust (or sphericity) axis. Note 
that the gluon jets, coming from hard gluon radiation, 
will most probably satisfy this condition. 

In this approach the basic questions are: How to iden- 
tify the jets? How to partition the phase space further 
after isolating out the individual jets? Let us try to an- 
swer these questions. 

The identification of jets from the final-state multipar- 
ticle system can be done using some algorithm, e.g. the 



FIG. 2: (o)-(c) The same event as in Fig. 1(6) plotted in 
momentum space, in z-y (a), x-y (6) and z-x (c) planes, re- 
spectively, (d) The number of jets in this event determined 
by different values of d cu t- 



JADE or the DURHAM || ones. In these schemes a 
distance measure is defined for counting the separation 
between two jets (or two particles). For example, in the 
DURHAM algorithm this measure is defined as 



V 1 j; (l-cosfl 



(4) 



which is related to the relative transverse momentum fct 
as 



(fct)* 



(5) 



In Eq.(4) Ei and Ej are the energies of the two jets (or 
two particles), 9ij is the angle between them and s is the 
square of cm. energy of the event. Note that we have 
used the notation dij instead of yij , currently used in the 
literature, to avoid confusion with the rapidity cut. Jets 
or particles with dij < d cu t are combined into a single 
jet. This procedure is repeated until all pairs of jets i 
and j satisfy dij > <i C ut- 

In principle, d cn t could be chosen arbitrarily. The re- 
sulting "jets" are, however, very different for different 
values of d cu t. For our present purpose, we want to cut 
the phase-space according to jets, so the "jets" should 
be directly observable, cf. Fig's. 1 and 2. Such directly 
observable jets will be referred to as "visible jets" . 

It has been shown in Ref. JTo| that the scale for visible 
jets is 5 < fc t < 10 GeV. For the case of */a = 91.2 
GeV, which we are taking as example, this range of k t 
corresponds to 10 -2 5 < d cu t < 10~ 19 . 

In order to be more apparent, we plot in Fig's. 2 (a)-(c) 
the view in momentum space of the same event, shown in 
Fig. 1(6) as lego-plot. The number of jets iVj et versus log 
d cu t in this event is given in Fig. 2(d). It can be seen from 
the figure that if we were using a value of d cut as small 
as 10~ 4 (or 10~ 3 ) then we would conclude that this is a 
10-jet (or 5-jet) event; while if we used a value of cZ cu t 



FIG. 3: The rapidity distributions (histograms) of 2-jet 
events determined by DURHAM algorithm with different values 
of dcut- Full triangles — that of the events with the bigger jets 
at right. Open circles — that of the events with the bigger 
jets at left. 



greater than 10 _1 then we would take this event as a 
2-jet one. Although all these values of d cut are equally 
possible in principle, it is very unlikely that there are 10 
(or 5) jets or only 2 jets in this event, cf. Fig's. 1 and 2. 
For our present purpose of constructing extended phase 
space through plotting sub- lego plots of individual jets, 
this event should certainly be considered as a 3-jet event. 
It can be seen from Fig. 2(d) that a value of d cut about 
equal to or a little smaller than 10~ 2 will combine the 
wrongly identified five jets into three. This is the d cu t 
value corresponding to visible jets. 

A further check of this statement could be obtained 
from the rapidity distributions of 2-jet events for four 
different values of d cu t, shown in Fig. 3 as histograms. In 
the same figures, the full triangles (open circles) are the 
distributions for those events in which the bigger jet is at 
the right (left). It can be seen from the figures that for log 
d C ut = —2.5 ~ —2.0 the full triangles manifest a shoulder 
at about y = 1.8, showing that at this value of d cu t a third 
jet starts to appear inside the big jet. This confirms the 
assertion that log d cut = —2.5 ~ —1.9 is the appropriate 
distance-measure for identifying the visible jets in e + e _ 
collisions at yfs = 91.2 GeV. For convenience, we will in 



and use a value 



the following refer to this measure as d£ut 
= 10-". 

Now let us turn to the second question: How to parti- 
tion the phase space further inside the identified jets? 

The first step is to take out the identified jets and 
draw sub-lego plots for them, forming an extended phase 
space together with the mother lego plot. In these sub- 
lego plots the jet axes are taken as the center to define 
new variables 



ff = 77 - JJjet, If' 



f - <Pjet- 



(6) 



As example, the sub-lego plots of the two identified jets 
in Fig. 1(6) are shown in Fig. 4. 

The pair of variables (J7j e t, fjet) determine the direc- 
tion of jet axis. A circular symmetry is expected to exist 
around this point in the rj-tp plane. Therefore, a natural 
variable for the description of the phase space region in- 



FIG. 4: 
Figl(fe). 



The sub-lego plots of the two identified jets in 




FIG. 5: (a) The distribution of r and (b) that of the corre- 
sponding cumulant variable r c . 



side a jet is the distance r in rj-tp plane from the particle 
inside the jet to the jet axis (77j e tj fjet), cf. Eq.(2). 

The distribution of r is shown in Fig. 5(a). It has a 
pick at about r — 0.3 and tends to zero when r increases. 
Choosing a certain value of Rq, e.g. Rq = 1-5, shown 
as solid circle in Fig's. 4, a phase space region can be 
defined as < r < Rq . A partition of this region using 
different values of r results in a series of rings, cf. the 
dashed circles in Fig's. 4. The distance Sr between the 
neighboring rings is the scale of the partition. 

It is interesting to check whether there is a scale in- 
dependence or anomalous scaling in the 1-dimensional r 
space, or in other words, whether the r-distribution is a 
self-similar fractal. This can be done using the standard 
technique [ pd| . 

The first step is to change the variable r into the cor- 
responding cumulant variable r c defined as fll2fl 



p(r)dr 




p(r)dr 



(7) 



The 



to eliminate the influence of nonflat distribution, 
resulting distribution of r c is shown in Fig. 5 (b). 

The second step is to divide the region < r < Rq (0 < 
r c < 1) into M bins and define the gth order normalized 
factorial moments nfm as [13l 



M 



(n m (n m -!)■■■ (n m - q + 1)) 



(8) 



where n rn is the number of particles in the mth bin, (■ • •) 
denotes the average over the sample of identified jets. 
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FIG. 6: The log F q versus log M plot. 

The advantage of NFM is that it can eliminate the statis- 
tical fluctuations provided the latter is of the Poinsonian 
type The anomalous scaling of nfm 

F q (M) oc M^ q M->oo (9) 

as the increasing of the partition number M or as the de- 
creasing of the phase space scale Sr = Rq/M, will signal 
the existence of fractal property — self-similarity. 



In Fig. 6 is shown the log F q versus log M plot in the 
region < r < tt |Q for q = 2,3,4,5. It can be seen 
that very good scaling is obtained for M — 2 - 40. This 
astonishing phenomenon is closely related to the color co- 
herence — angular ordering in parton fragmentation and 
is worthwhile further investigation, both experimentally 
and theoretically. 

In this letter the proper way for the phase space par- 
tition in high energy e + e _ collisions is analysed in some 
detail. It is shown that the partition of phase space ac- 
cording to the value of rapidity, popularly used in hadron- 
hadron collisions, is inappropriate for the study of e + e~ 
collisions. The proper way in the latter case is to iden- 
tify the visible jets and take them as objects for detailed 
study, forming an extended phase space. The value of 
the distance-parameter d cut of the DURHAM algorithm, 
which is able to identify the visible jets, is discussed in 
detail and is found to be about d^t = 10~ 2 ' 2 (at cm. 
energy y/s = 91.2 GeV). 

A new variable r is introduced to further partition 
the phase space inside jets, which possesses a very 
good anomalous scaling property, showing that the re- 
distribution inside jets is a self-similar fractal. This spec- 
tacular phenomenon, being connected with the underly- 
ing dynamics, is worthwhile further theoretical and ex- 
perimental investigations. 
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